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Abstract. Within the quantum mechanical treatment of the decay problem one finds that at 
late times t the survival probability of an unstable state cannot have the form of an exponentially 
decreasing function of time t but it has an inverse power-like form. This is a general property 
of unstable states following from basic principles of quantum theory. The consequence of this 
property is that in the case of false vacuum states the cosmological constant becomes dependent 
on time: A — Abare = A(t) — Abare ~ 1/t 2 . We construct the cosmological model with decaying 
vacuum energy density and matter for solving the cosmological constant problem and the 
coincidence problem. We show the equivalence of the proposed decaying false vacuum cosmology 
with the A (t) cosmologies (the A(t)CDM models). The cosmological implications of the model of 
decaying vacuum energy (dark energy) are discussed. We constrain the parameters of the model 
with decaying vacuum using astronomical data. For this aim we use the observation of distant 
supernovae of type la, measurements of H(z), BAO, CMB and others. The model analyzed is 
in good agreement with observation data and explain a small value of the cosmological constant 
today. 


1. Introduction 

The nature and origin of an accelerating expansion of the Universe is the basic problem of modern 
cosmology. The most natural explanation of the acceleration of the Universe seems to be the 
cosmological constant parameter interpreted as vacuum energy. Although such an explanation is 
simple it brings a conundrum that the values of the cosmological constant required by quantum 
theory (~ 10' 1 GeV 4 ) and obtained from type la supernovae observations (p\ = ~ 10 -47 

GeV 4 ) differ about 100 orders of magnitude [T]. This problem is called the cosmological problem. 

One approach toward the solving this problem is to consider the cosmological model with the 
time varying cosmological constant A(t), where t is the cosmological time. The parametrization 
for decaying vacuum is usually taken by hand to the cosmological model [21 [31215]. Our idea 
is to derive the parametrization of decaying vacuum energy directly from the first principle, 
namely from the quantum mechanics, then to construct the cosmological model and to test it 
statistically using astronomical data. Therefore our approach to decaying vacuum cosmology is 
not purely phenomenological and is motivated by the fundamental theory of quantum mechanics. 

1 Talk given at Seventh International Workshop DICE 2014: Spacetime - Matter - Quantum Mechanics ... news 
on missing links, Castiglioncello (Tuscany, Italy), September 15-19, 2014 and submitted to the Proceedings of 
this conference. 


Decaying false vacuum states from the point of view of the quantum theory of unstable 
states evolve in time. From basic principles of quantum theory it is known that the amplitude 
A(t), and thus the decay law VM(t) = |A(f)| 2 of the unstable state | M), are completely 
determined by the density of the energy distribution function u>(E) for the system in this 
state: A(t) = J^° . u(E) exp [—^Et\dE, where u(E) > 0 and ui(E) = 0 for E < E m i n . 
From this last condition and from the Paley-Wiener theorem it follows that there must be 
\A(t)\ > A\ exp [—A 2 t q ], for \t\ — > oo. Here A\ > 0, A 2 > 0 and 0 < q < 1. This means 
that the decay law Vm( t) of unstable states decaying in the vacuum can not be described by an 
exponential function of time t if time t is suitably long, t —> 00 , and that for these lengths of time 
Vu{t) tends to zero as t — > 00 more slowly than any exponential function of t. It appears that 
these deviations from the exponential decay law at long times affect the energy of the unstable 
state and its decay rate at this time region and thus they affect the energy of the unstable false 
vacuum states at these times. It is shown in 013 that at transition times t ~ T, where T 
denotes time when contributions of the exponential part of the survival probability and of its 
late time non-exponential part are equal, the instantaneous energy of the false vacuum states 
fluctuates and at late times, much latter then transition times, t T, it tends to the energy of 
the true vacuum state as 1/f 2 for t —> 00 . The asymptotically late time behavior of the energy 
of the system in the false vacuum state is given by the following relation 

2 

£f se (f) ~ E$ ue ± ^ • • • , for t » T. (1) 

This means that in the case of such false vacuum states the cosmological constant becomes 
time depending. The standard relation is pQ me = Pbare = ■ So the fluctuations of 

PQ alse (t) = EQ alse (t)/V at the transition times region and the asymptotically late behaviour 
of p() alse (t) at t T mean that identical behaviour of A have to be observed at these times. 

There are two possible scenarios. In the cosmological terminology if the universe is in a false 
vacuum state then the cosmological constant is time-dependent in the following way 

3/3 

A(t) = A bare + -p, for t » T. (2) 

In the model with decaying dark energy and matter, as in most cosmological models, we act on 
assumptions that the Universe is homogeneous and isotropic at the large scale (the cosmological 
principle), the structure and evolution of the Universe is governed by the Einstein theory of 
relativity, and the source of gravity is described by energy-momentum tensor for matter and 
decaying vacuum. It is also assumed for simplicity that the Universe is flat. 

Then cosmological evolution is determined by the scale factor a as a function of the universal 
cosmological time t, which satisfies the Friedmann equation 
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where matter density fulfils the conservation condition 

T a0 .p = + A(t)g a0 = 0 . 

This condition for the flat homogeneous and isotropic Universe has the form 


( 3 ) 

( 4 ) 


Pm T 3 Hp m — A(f) 


( 5 ) 


In the case of the ACDM model (/3 
after substitution a 3 = x 2 


: 0) it would be useful to rewrite ^ to the new form 

3 3 . 2 

^Pm,0 T ^ A b are® ■ (6) 



In cosmology we use instead of density of fluids dimensionless parameters called density 
parameters 

,U 3 Hi 

where pi is energy density of the fluid i, Hq is the Hubble function H = a/a and the index 0 
denotes the present epoch. 

Therefore eq. ([ 6 ]) assumes the form 


dx\ 2 9 9 2 

JJ ) ~ 4^ m ’° + l^bare.O® 


( 7 ) 


where H A , = 


bare 3 Hf, 


and r = \Ho\t. 


Then we obtain the solution in the form 
x(t) = 


O \ >/2 

“m,0 


Ha 


bare^O 


Sillh H 0 t \ . 


After dividing both sides of eq. @ by 3 Hj we obtain the relation 


( 8 ) 


) — ^m, 0 a 3 + ^A b 


(9) 


This relation can be rewritten in the terms of redshift z &s 1 + z = a 


— 1 


(rj- ) - H mj0 (l + z) 6 + H Abarei0 . 


( 10 ) 


Today (z = 0 and H = Hq) density parameters satisfy the constraint relation fi A ,o = 1 — H m ,o- 
This relation is fundamental for cosmography which analyzes the observational effects of photons 
propagation along the zero geodesics. 

For the model with vacuum decaying 
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j ~ ^ m ( Z ) + ^Abare + T(z) 2 


( 11 ) 


where 


T( Z )=- r 

J oo 


dz 


(1 + z)H(z) 


is the age of the Universe up to the redshift z. 

For estimation of the model parameters we assume that 


H{z) = Ffo(l + z) 7 , 7 = const > 0 


which gives 
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T(z) = -^-(1 + z )" 7 oc H- 1 = t H , 


iHo 


where tn is a Hubble scale. Then 


+ HA bare + ^"77 2 - 


( 12 ) 







The acceleration equation has the form 


H = - 6H 2 , 


s = -(i~P). 


After changing the variable Mawe obtain the acceleration equation in the form 

dH do A bare 


da dt 


-SH¬ 


OT 


dH 1 / A bare o 
lO = Wa [ ^ ~ 6H 


The first integral of eq. (15) is 


H 2 (a) =[H 2 - 


■A-bare ® 
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~ + 


-A-bare 


(13) 

(14) 

(15) 

(16) 
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From this relation we obtain H 2 (z) formula which will be used in model parameter estimation 


(10 ~h Qmil+z)2S+ h QAbm 

where (further let write down Abare = A) 

2 

ilm.O + f^A,0 = — 5 - 

When 6 = 3/2 the model reduces to the standard ACDM model. 
Finally, we obtain the following model 


(17) 


(~ 

\Ho 


^ m >° (x + z j 3 (fl m , 0 +fi A ,o) + 

^m,0 + ^A,0 ^m,0 + ^A,0 


(18) 


After the generalization of the corresponding substitution 0 3 (^m,o+s*U,o) = x 2 ) we obtain the 
expression for the scale factor in following form 


a{t) = 


^111,0^ 3 ( n m,0 + n A,0) 

^A ,0 


sinh ( -^^AflH 0 t 


3 (. n m ,0 + n A, 0 ) 


2. Data 

For the estimation of parameters of the model (18) we used the observations of SNIa data, BAO, 
CMB, measurements of H(z), Alcock-Paczyhski test. 

First, we consider the SNIa data. The likelihood function is 


In TsNIa 



(19) 


where the summing is over the SNIa sample; the distance modulus / u obs = m — M (where m is the 
apparent magnitude and M is the absolute magnitude of SNIa stars) and /u th = 51og 10 Dl + 25 
(where the luminosity distance is Dl = c(l + z) J 0 ~ jjfzrj and a is the uncertainties. The Union 
2.1 sample of 580 supernovae was used as the data here [8]. 



















The BAO (baryon acoustic oscillation) data were taken from the Sloan Digital Sky Survey 
Release 7 (SDSS R7) dataset which consists of 893 319 galaxies [5]. The likelihood function is 
given by 


In Tbao 


r s (zd) 

D v (z) 


~d(z)^ 




where r s (zd) is the sound horizon at the drag epoch and z = 0.275, d(z ) = 0.1390, a = 0.0037 

[El- 

Planck observations of cosmic microwave background (CMB) radiation were also used [11] . 
With addition information on lensing from the Planck and low-£ polarization from the WMAP 
(WP), we obtain the combined likelihood function in the form 


InLcMB+lensing+WP = ~ - X° bs ), (20) 

ij 


where C is the covariance matrix with the errors, x is a vector of the acoustic scale l a, the shift 
parameter R and f l^h 2 where 


Ia 



dz' 

Wz 0 


R = 



dz' 


( 21 ) 

( 22 ) 


where z* is the recombination redshift. 

The Alcock-Paczynski test is the comparison of the radial and tangential size of an object, 
which is isotropic in the correct choice of model mm • The likelihood function is independent 
on the Hq parameter and has the following form 


In Lap = 


1 (AP th ( Zi ) - AP ohs (zi)Y 


o 


(23) 


where AP(z) th = / Q ~ jffY) an ^ AP(zi) ohs are observational data [14, [TSJ 16, T7, 18,, 19, 20, 

mm- 

At the end it is also valuable to add the constraints on the Hubble parameter, i.e. H(z = 
0) = Hq. Data of H(z) for samples of different galaxies were also used [231 EH, 25) 


In L H{g) = -\YY 


i= 1 


H(zi) obs - H{zif 


(Ti 


The final likelihood function for the observational Hubble function is 
Ltot = T'SNIa-T'BAO -^CMB+lensing+WP -^AP ^H(z) ■ 


(24) 


(25) 


3. Estimation of the model 

To estimate the model parameters we use our own code CosmoDarkBox implementing the 
Metropolis-Hastings algorithm [26, [27J and using the Pade approximants |28j for the calculation 
of the likelihood function mm- 

We use observation data of 580 supernovae of type Ia, selected subsets of the data points 
of Hubble function, the measurements of BAO from SDSS+2dSGRS. We also use data for the 












Table 1. The values of estimated parameters - best fitted values of parameters with errors. 
We consider two cases for errors estimation, first Hq is assumed as 68.22 km/(s Mpc), then H m 
is assumed as the best fit value. 





parameter 

best fit 

68 % CL 

95% CL 

Case 1 

- assumed Ho 

= 68.22 

1 

r^m 

0.0068 

0.2926 

+0.0261 

-0.0276 

+0.0251 

-0.0239 

+0.0419 

-0.0461 

+0.0416 

-0.0388 

Case 2 

- assumed H m 

= 0.2926 

1 H m — Ha 
Hq 

0.0068 

68.22 

+0.0200 

-0.0203 

+0.92 

-0.91 

+0.0324 

-0.0331 

+1.50 

-1.45 



Figure 1. The likelihood function of two model parameters (1 — H m — Ha, H m ) with the marked 
68 % and 95% confidence levels. The value of Hubble constant is estimated from the data as 
best fit value Hq = 68.22 km/(s Mpc) and then the figure is obtained for this value. 


application of the Alcock-Paczynski test - 18 observational points. At last, we estimated model 
parameters with CMB data from Planck, low-1 polarization from WMAP and lensing from 
Planck. 

The results of statistical analysis are represented in Figures 0B The values of estimated 
parameter are completed in Table 1. In Figures [l] and [2] it is shown the likelihood function 
with 68% and 95% confidence levels projection on the (1 — fi m — HA,H m ) plane and the 
(1 — fi m — Ha, Ho) plane, respectively. In figures [3] and [4] the intersections with respect to 
fixed Hq and Ha for H m : and fixed Hq and fi m for Ha have been presented, respectively. 

We obtain, that the effect of decaying vacuum is rather small and errors one order higher. 
Therefore the effect cannot be detected from available observations data. Second, the errors for 
parameter 1 — H m — Ha are the same order like fi m , Ha and we are looking for small parameter 
close to zero. 

4. Conclusion 

(i) The cosmological model with decaying vacuum explains while the value of cosmological 
constant is so small — possible solution of the cosmological constant problem. 

(ii) The model with decaying vacuum (cosmological constant) was tested using astronomical 








Figure 2. The likelihood function of two model parameters (1 — D m — Q\,Hq) with the 
marked 68% and 95% confidence levels. The value of D m is estimated from the data as best fit 
value 0.2926 and then figure is obtained for this value. Note that the effect of decaying vacuum 
appears as the value of the parameter 1 — D m — Da is different from zero. 
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Figure 3. Diagram of PDF for parameter 1 — D m — Da obtained as an intersection of the 
likelihood function. Two planes of intersection likelihood function are Hq = 68.22 km/(s Mpc) 
and D m = 0.2926. The planes of intersection are constructed from the best fitting value of the 
model parameters. The maximum of PDF is reached for 1 — D m — Da = 0.0068. 


data (SNIa, BAO, CMB, H(z)). 

(iii) The evidence of decaying dynamical vacuum effect for the current Universe is equivalent to 
the sum of D m] o + Da,o 7 ^ 1- For the standard cosmological model this sum is equal one. 

(iv) The value of sum of D m and Da is close to 1 (the obtained value is 0.993) — the effect of 
decaying vacuum is very weak. 

(v) The cosmological models with decaying vacuum be can treated as an extension of the 
standard cosmological model. This model is in a good agreement with astronomical data, 
and it offers the solution of the cosmological constant conundrum. 
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Figure 4. Diagram of PDF for parameter D m obtained as an intersection of the likelihood 
function. Two planes of intersection likelihood function are Hq = 68.22 km/(s Mpc) and 
1 — Q m — Da = 0.0068. The planes of intersection are constructed from the best fitting value of 
the model parameters. The maximum of PDF is reached for D m = 0.2926. 
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